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1 Introduction 



Over the last decade, following the paper of Bartels and Loewe[l], diffractive leptoproduc- 
tion of vector mesons has been investigated within the framework of perturbative QCD 
(pQCD). In particular, it has been shown[l,2] that the DIS cross sections can be calculated 
in pQCD provided that both the energy and photon virtuality are large enough. These 
cross sections are proportional to {xG{Q'^ , where xG{Q^, x) is the gluon distribution 
within the nucleon target. Futhermore, Donnachie and Landshoff[3] have shown in a non 
perturbative Pomeron approach, that both the initial photon and the produced vector 
meson have to be longitudinally polarized to realize the leading behaviour of the cross 
section i.e. ^ oc This result has also been confirmed in the BFKL approach to the 
Pomeron structure in pQCD[4]. A new insight into the problem has been suggested by 
Ryskin[5], who proved that J/ip leptoproduction can be calculated reliably in the leading 
log approximation (LLA) of pQCD, and that the non-perturbative effects coming from 
large distances, can be factored out in terms of the wave function of the produced vector 
meson at the origin. In addition, Ryskin pointed out that the pQCD calculation is also 
valid for leptoproduction of light vector mesons. For heavy vector mesons the calculation 
can be safely continued to — 0. An important contribution in the construction of this 
formalism appears in the paper of Kopeliovich et.al.[6], where the wave function for the 
light vector mesons is evaluated using the constituent quark model. Finally, we mention 
the contribution of Brodsky et.al.[7] who showed how to factor out the long distance ef- 
fects for the case of light vector mesons in terms of the light-cone qq wave function of the 
produced vector meson [8]. 

The goal of this paper is to study the shadowing (screening) corrections (SC) for lepto 
and photoproduction of vector meson. To this end we shall generalize the formalism 
developed by Mueller [9] for the gluon and quark densities. Technically, we shall utihze 
the idea of a r±- representation for the SC suggested in Ref.[10] and further discussed in 
Refs.[6,9,ll,12]. The quantity we wish to calculate is the damping factor defined in our 
previous publication [13] as: 



where t = —qj_ denotes the square of momentum transfer in the reaction 7*p — > Vp . 

The structure of our paper is as follows: In section 2 we define our notation and 
numerical coefficients pertinent to the DIS production of vector mesons without includ- 
ing SC. Basically, we reproduce the well known results of Refs.[5,7], doing calculations 




inth SC] 



(1) 



dzh:£fm[ without SC] 



t=0 



1 



in the r_L-representation ( r± denotes the transverse separation between the quark and 
antiquark) . 

We will show that the results hold, not only in the leading log of the GLAP evolution 
equations [14], but also in the leading In- approximation of pQCD, as was pointed out 
by Ryskin[5]. In section 3 we extend the formalism suggested by Levin and Ryskin[10] 
and Mueller [9] to the case of a qq pair with a definite value of r± transversing the target 
(be it a proton or nucleus). We will show that for vector meson production we obtain 
a closed expression for the damping factor. Our result does not depend on the large 
transverse distances, allowing us to calculate the SC to a good theoretical accuracy in 
pQCD. In the fourth section we prove that the damping factor can be rewritten using 
the experimentally measured value of To this end, we use the leading log ^ 

approximation of pQCD with a variety of shadowing corrections. In section 5 we discuss 
the SC for the gluon distribution, and investigate the uncertainties due to large distance 
contributions which can be important for this case. Section 6 contains the numerical 
evaluations of our formalism, which arc compared with those obtained for the case with 
no SC. Our results are then compared with the experimental data, mostly from HERA. 
A summary and final discussion are presented in section 7. 



2 Vector meson production in DIS without shadow- 
ing corrections 

2.1 Notations 

In this paper we will use the following notations (see Fig.l): 

denotes the virtuality of the photon in DIS, my is the mass of the produced vector 
meson and m is the target mass assumed to be a proton unless specified otherwise. 

X — ^_±!!V^ where ^/s is the cm. energy of the incoming photon- proton system. 

0^= z(l — z) Q"^ + rriQ, where z is the fraction of of the photon energy carried by the 
quark, mq is the mass of the current quark. 

— * 

k± denotes the transverse momentum of the quark, and f± the transverse separation 
between the quark and antiquark. 

— * 

b± is the impact parameter of our reaction which is the variable conjugate to the 
momentum transfer {q±). 
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li± denote the transverse momentum of the gluons attached to the quark- antiquark 
pair (see Fig.l). 

We will use the evolution equations for the parton densities in the moment space. For 
any function f{x), we define a moment f{u!) 

fiuj) = f' dx X- f\x) . (2) 
Jo 

Note that, the moment variable uj, is chosen such that the = moment measures the 
number of partons, and the u = 1 moment measures their momentum. An alternative 
moment variable N defined such that N = + 1 , is often found in the literature. 

The X - distribution can be reconstructed by considering the inverse Mellin transform. 
For example, for the gluon distribution it reads 

xG{Q\x) - J^d^ 9{Q',^) ■ (3) 

The contour of integration C is taken to the right of all singularities. The solutions to 
both the GLAP and the BFKL equations have the general form 

giQ^u) = ^(a;)e^H^"«\ (4) 

where 7(0;) denotes the anomalous dimension, which in the leading In^ approximation of 
pQCD, is a function of ^ and can be presented as the following series[15] 

^ ^ + 0(4) . (5) 

Our amplitude is normalized such that 

(Jtot = 4:71 lmf{s,0) . 
The scattering amplitude in 6_L-space is defined as 

a(sM) = ^ / d'q^e-^^^--'- f{s,t^-ql) . 
In this representation 

atot = 2 d%_L Im a(s, b±) , 
aei = j d%± I a{s,b±)f . 

In general we have attempted to adjust our notation and normalization to agree with 
those of Brodsky et.al.[7]. 
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2.2 The amplitude in the r± representation 



This approach was originally formulated in Ref.[10] and has been carefully developed in 
Ref. [9]. While the boson projectile is transversing the target, the distance r±_ between the 
quark and antiquark can vary by an amount Ar± oc i? where E denotes the energy 
of the pair in the target rest frame, and R is the size of the target (see Fig.l). The quark 
transverse momentum is A; i oc — . Therefore 



Ar^ oc < , (6) 

which is valid if 

■ s > 2m i? , (7) 
where s ~ 2mE. The above condition has the following form in terms of x 

X < ^— . (8) 
2mR ^ ' 

Hence, at small values of the transverse distance between the quark and antiquark is a 
good degree of freedom [9, 10, 16], and the interaction of a virtual photon with the target 
can be written in the form ( we use the notation of Ref. [7]) 

Mf ^ ^fl^f £ J ^ ^'\Q',r^,z) a(rx,gi) [^''{r^,z)]* , (9) 

where denotes the polarization of the quarks. Note that both 7* and V are longitu- 
dinally polarized. denotes the wave function of the produced vector meson, and we 
anticipate that the value of rj_ that dominates the integral in Eq.(9) will be small (of the 
order of oc ^). The form of \E''^*, the wave function of the longitudinally polarized 
photon, has been given in Ref. [17], i.e. 

*^*(r±,^) = Qz{l-z)Ko{ar±) . (10) 

a{r±, qY) is the cross section of the qq pair with a transverse separation r± which scatter 
off the target with momentum transfer q±. We first evaluate it at q± — 0. 



2.3 a{r±,ql) at q± = 

The expression for a{r±, qjj at q± — 0, was first written down in Ref. [10] (see Eq.(8) of 
this paper) . It turns out that a can be expressed through the unintegrated parton density 
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(j) first introduced in the BFKL papers [18] and widely used in Ref.[2]. The relation of this 
function to the Feynman diagrams and the gluon distribution can be calculated using the 
following equation: 

as(g') • xG{Q\x) = f dll asill) </.(/!, x) . (11) 
Using the above equation we reproduce the results of Ref.[10], which reads 

where = ^'^aiy ''^^ and Cp = ^^^T' ^® evaluate this integral using Eqs.(3,4) and 
integrate over the azimuthal angle. Introducing a new variable ^ = r^Z^, the integral can 
be written in the form 

^(-^' - ^ Jo ^^^^ ^^^^ io ^ ^ (0^ • ^''^ 

Evaluating the integral over ( (see Ref.[19] 11.4.18) we have 

In the double log approximation of pQCD, where 7(0;) <^ 1, the cross section for Nc — 
Nf — 3 reads 

-(r.,?i) = ^ rl ( xG«-'^-(A j . (15) 

This result coincides with the value of the cross section given in Refs.[6,20], (if we neglect 
the factor 4 in the argument of the gluon density). We checked that Eq.(2.16) of Ref. [7] 
also leads to the same answer, unhke the value for a quoted in Ref.[7] (see Eq.(2.20)) 
which differs from ours by a factor of 2. 

We do not want to make use of the double log approximation, which corresponds to the 
first term in the expansion given in Eq.(5), and prefer to use the complete series for 7(0;). 
Indeed, all diagrams in which an extra gluon goes from the bottom to the top as well as 
from left to right of the diagram of Fig.l (see also Fig. 2) do not produce an extra power 
of In -. They only give corrections of the order as either to the anomalous dimension, or 
to the coefficient function. Calculating in the In^ approximation, considering ^ 3> as 
in the region of small x, we can safely use the value of the anomalous dimension given in 
Eq.(5). From Eq.(14) one can directly obtain the simple answer in another limiting case. 
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namely, when 7(0;) — > |. This hmit corresponds to the solution of the BFKL evolution 
equation[18] in the diffusion region when In 4- < \/ttsln-. Neglecting the deviation of 
7 from I we get 

a(r„,i) = ^vrVi (xG^ (A,) j . (I6) 

In the intermediate range of and x, we cannot obtain a simple formula for a in terms 
of the gluon density. We stress that a difference of a factor of two in a does not mean 
that the cross section of vector meson production will change by a factor four. The BFKL 
density behaves as r±F{lnr'j_, x) , and the integral over r± will be different for the double 
log approximation and for BFKL. We will consider these two cases and will integrate over 
r_L to see how big the difference is in the final answer. 



2.4 The amplitude of vector meson production at q± = 

We wish to calculate the amplitude for vector meson production, using Eq.(9). To do this 
we use the integral representation for the McDonald function [19] and rewrite Eq.(lO) in 
the form 

Substituting Eqs.(14,16) in Eq.(9) we obtain 

(iV_L Qz{l-z) dt ( 1 o?'r\ \ 



t 



ic27^^ ''"'M' {T{2-'y{uj)y 

Anticipating that the typical value of r± in the integral will be small {r± oc ^), we 
can replace the hadron wave function by its value at r±— 0. Introducing a new variable 

2 2 

C = we integrate over r^^ and t. It is easy to see that the integral corresponds to 

C ~ 1 and t ~ 2. This means that the typical value of r'j_ is of the order of 4^ ^ 1- This 
justifies our assumption that the wave function of a produced hadron enters at small r±. 
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The final general answer is 



Mf-J^fT.f'dz ^^ii^ ^ I ^ g{u) (19) 



Ai A2 



1 fi 

- (a2)^(-) r(l + 7M) r(l - 7(cu)) [^^(rx = 0, z)]^ 



In the case of the GLAP approach we can simplify Eq.(19) by considering 7(cj) ^ 1, 
while for BFKL the simplification occurs when 7(0;) —>■ |. One can see that there is a 
difference of 7r/2 in the amplitude, but the final estimate is obtained only after integration 
over z. 

We can approximate the integration over z by putting z = ^ in the Ina^ 
dependence of the density. This gives 

Mj = CQ as (xG«^-^^(a^(^ = i),a;)) £ dz ^il-l^[^^(0, ;.)]* , (20) 

for the GLAP equations, while for the BFKL we get 

Mf = CQas \ (xF^'^\\na'{z^\),x)) £ dz ^^^[*^(0, .)]* , (21) 

where we have replaced the BFKL gluon distribution by xG^^^^{a'^, x) = ^ • F(lnr^, x) 
as discussed above. In the above equations we have replaced all constants by the coefficient 
C. 

The function ^'^(O, ^) is just the function which has been introduced in Ref.[7] (see 
Eq.(2.25) there). Using the asymptotic form of the z-dependence, namely oc z{l — z) 
for the case of massless quarks, we have: 
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= [xG^'^^i^,x)j , (22) 

= (^„x)) . (23) 

Consequently, the difference in the value of the amplitude turns out to be rather small, 
about 25%. We note that the BFKL equation gives the anomalous dimension close to |, 
which only generates a ^ behaviour of the amplitude, while the GLAP equations give 
^ for the same value. In the nonrelativistic case the difference is 7r/2 in favour of the 
BFKL equation. 
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2.5 The b± dependence of the amplitude 



To deal with the SC we need to know the ampHtude not only at q±=0, but also at all values 
of momentum transfer, so that we can calculate the profile function of the amplitude in 
impact parameter space. The gluon density has a week dependence on q± at small values 
of q±, both in the double log approximation (see Ref.[2]) as well as for the BFKL approach 
(see Ref.[20]). Therefore, the gfj_-dependences come from the form factor of the qq pair 
with a transverse separation r±, and the form factor of the proton target. 

The target form factor is not treated theoretically in pQCD, and that for our purpose 
we assume the exponential parameterization for it will suffice, namely 

FMl) = e"^ ■ 

If we put the Pomeron slope a'=0 the slope B can be extracted from the experimental 
data on hadron-hadron collisions. Namely, B = B^[a' = 0), where B^i is the slope of the 
differential cross section of proton-proton collision. For the numerical estimates we use 
phenomenological information to extract the value of B [21]. The value of B defined in 
this way, is very close to the one obtained from the proton electromagnetic form factor. 

The form factor of the qq pair with a transverse separation r± is equal to 

where ki (/cQ denotes the momentum of the i quark before and after the collision. Each 
of the wave functions has an exponential form, and a simple sum of different attachments 
of gluon lines to quark lines gives 

F,M) = e'^{l - e^^---}. (25) 

We have absorbed the last factor in the expression for the cross section, while the first one 
gives the q± dependence of the qq form factor, which after integration over the azimuthal 
angle has the form 

F,M) = •^o(^) . (26) 
To calculate further, we need the profile function in b± space, which is defined as 

"^^^'^ ^ i / ^'^^ ^^^^'^ ^^^^^'^ ■ ^^^^ 
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An exact computation gives 



To simplify the calculation we replace the above function by 



S(bl) = ^ e-* , (29) 



where 

B' ^ B ( I + ^] B ( 1 + ^ 



AB J V ■ a^B 

Eq.(29), in the small b\ expansion, has the same radius as in Eq.(28). 

3 Shadowing corrections to diffractive 
leptoproduction of vector mesons 

3.1 Percolation of a qq-pair through the target. 

To calculate the shadowing corrections wc follow the procedure suggested in Refs.[9,10,13,21]. 
Namely, we replace a{r±, g^. = 0) in Eq.(9) by 

a^^(rx) = 2 Jd%^{l - e-^-('^-'^i=o)^(^i)) . (30) 

The above formula is a solution of the s-channel unitarity relation 

2Ima{s,b±) = |a(s,6x)|' + Gin{s,b±) , (31) 

where a denotes the elastic amplitude for the qq pair with a transverse separation r±, and 
Gin is the contribution of all the inelastic processes. The inelastic cross section is equal 
to 

a^n = J d'b^GUsM) = jd%^ (l - e- '^('•-^i=°) ^(^i) ) . (32) 

We assume that the form of the final state is a uniform parton distribution that follows 
from the QCD evolution equations. Note that we neglect the contribution of all diffraction 
dissociation processes to the inelastic final state (in particular to the "fan" diagrams which 
give an important contribution[2]), as well as diffraction dissociation in the region of small 
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masses[21], which cannot be presented as a decomposition of the qq wave function. We 
evaluate this input hypothesis in the next section. 

In the language of Feynman diagrams, Eq.(32) sums all diagrams of Fig.3 in which 

the qq pair rcscatters with the target, and exchanges "ladder" diagrams, each of which 
corresponds to the gluon density. This sum has already been performed by Mueller [9], 
and we will only comment on how one can obtain the result, without going into details. 

The simplest way is to consider the inelastic cross section (see Refs.[2,9]), which has a 
straightforward interpretation through the parton wave function of the hadron (see Fig.4), 
as all partons are produced on the mass shell. In leading In^ approximation we have two 
orderings in time (see Fig.4): 

1. The time of emission of each "ladder" by the fast qq pair, which should obey the 
obvious ordering for n produced "ladders" 

ti > t2 ... > ti > ti+i ... > tn . (33) 

2. Each additional "ladder" which should live for a shorter time than the previous one. 
This gives a second ordering 

ti-t[ > ^2-4 > - > ti-t'i > ... > tn-t'^. (34) 

Each "ladder" in the leading log approximation has the same functional dependence on 
ti — t'^, which we denote as a{t — t',r\). This fact allows us to carry out the integration 
over ti and ti —t'i, which gives for n emitted cascades 

^ j^f dt, P d{t^ - t[) a'\t^ - t[, rl) . (35) 

For the case of the qq pair, both last integrations are not logarithmic, and in the LLA we 
can safely replace the above integral with 

= a"(t-t',ri) , (36) 

t — t' is of the order of — due to the uncertainty principle, where q\\ = ^ ^"^v _ Qj^g 

is compensated by the number of possible diagrams, since the order of the "ladders" are 
not fixed. Therefore, the contribution of the nth "ladder" exchange gives 

<Tn ^ - (7"(x,ri) . (37) 

77/. 

Applying the AGK cutting rules [22], we reconstruct the total cross section which results 
in Eq.(32). 
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3.2 Damping factor for diffractive production of vector meson 



Substituting a'^'^' of Eq.(32) in Eq.(9) and using a in the form of Eq.(14) and the rep- 
resentation (17) for the wave function of the photon, we can easily estimate the general 
term of the expansion with respect to the power of a. It has the form 

^n_^ rdt_f 1,^ , a'rl,\ i-ir-' fasACFTT 



rn ^ f, - Z r°° dt / I , 



t 7 n\ \ J 

(38) 

^ f diJ, r(l + 7(^.)) r(l-7(a;,)) 



(r(2-7(a;.)) 



[*^(rx = 0,z)]* I dV^"(6l) . 
Taking the integrals over r^, i and hx_ we obtain 



- ^ ^ k '^'^1^^ \B\Nl - 1) ) ^^^^ 

du, . . r(l + 7(c^,)) r(l-7(a;,)) 



(r(2-7M)^ 



r\l + n-j:^{u;,)) (-) (a^)^^^'^') [*^(r^ = 0,^)]^ 



In the double log approximation of pQCD, 7(cjj) ^ 1 and we derive a very simple 
formula (neglecting 7). Taking the integral over uji we have 

/■! ^ g2:(l-z) (-l)-in! 



VS'(Ar2_i) a2 
Finally for Mf we have 



Jo 2a^ n ^ ^ 

^^"^ \asxG{a\x)\ [*^(rx = 0, ^)r . 



-CB'tl" dz^^^^ {-ir-\n - 1)! (41) 



n=l 



n 



CpT^ 4 2 \\ \iT,V, 



asxG{a\x)\ {r± = 0, z)] 
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The above series can be written in terms of the analytical function Ei (see Ref.[19]), 
namely 

-f ^CB' r dz9^^^ Ml)e'^ [^^{r^^O,z)]\ (42) 

JO Qj rv 



where (for Nc= 3) 



N = I • xG{a^, x) . (43) 



3 B' a? 

Using the above equation we have for the damping factor 

U',dzQ^K,^>y{r^^Q,z)Y 

For the case of heavy mesons we can use the nonrelativistic approach 2; — > |, and the 
above formula has the very simple form 

= ^^^^ . (45) 

Kg 

The behaviour of the damping factor can be easily found using the well known property 
of El. Namely, at ^ 1, we have 1 while at ^ 1 the damping factor vanishes 

as oc ^" 2'' ■ The behaviour of the damping factor as a function of k is given in Fig. 5. 



4 The relation between DIS difFractive production of 
vector mesons and F2(Q^,x). 

The expression for the F2 structure function was derived by Mueller [9] and in our notation 
has the following form 

E(Q^x) = I' dzj^ *r(Q^rx,^) <7^^(rx) [*r(g^rx,^)]* , (46) 

where 

= (47) 

(e^) denotes the average quark charge, which is equal to | (^) for three (four) active 
flavors respectively. 
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represents the wave function of the virtual photon with transverse polarization. 
, as needed for our calculations, has been given in Ref.[9] (see also Refs. [17,23]) i.e. 

*^(Q^rx,z) = aK^{ar^) . (48) 

However, as was shown in Ref.[9], within the LLA of pQCD, after integration over z in 
Eq.(46), we can safely replace this function by 4-- Finally Mueller's result reads (for 

iV, = Nf= 3) 

S Q^x = - / ^ . 49 

Comparing this equation with the general Eq.(9) we conclude that 

Using the notations of Ref. [7] and our normalization of the amplitude, we can rewrite the 
expression for the amplitude in the form 

/(7* ^V) = ^ f' dzQz{l-z) r rldr^ Koiar^) (51) 
V2 Jo 



ain^ / dz^^{z) ' 

where '{'^ {z) = \l/^(r_|_ = 0,z). 

The atractive feature of the above formula is that it takes into account all shadow- 
ing corrections in the LLA of pQCD, and expresses the cross section for vector meson 
production directly through the experimental observable. We note that without SC this 
formula is trivial and follows directly from the GLAP evolution equations in the LLA of 
pQCD (see, for example, Ref. [24]). The shortcoming of this formula is obvious: we cannot 
perform the integration over r± as has been done in Eq.(44). We think that this formula 
will be especially useful for reactions with a nuclear target, since it absorbs all SC which 
are essential in this case. 

Using Eq.(51) we can rewrite the expression for the damping factor in the form 



dzQzil-z) Jrldr^Mar^) , ,„ ^ ]f£^}' 



{J,'dzQz{l-z)Jrldr^ K,{ar^) f0l-^} 



2 



^We thank A. Mueller for a discussion with us on this problem. 
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In the case of a nuclear target, we have 



dln-kr dln-^ 

rj_ rj_ 

where F2{N) is the structure function for a DIS with a nucleon. The last remark makes 
the whole expression of Eq.(52) useful for a nuclear target, since all values that enter the 
formula can be measured experimentally. In the case of a nucleon target, we have to rely 
on the solution of the GLAP evolution equations to estimate the denominator. 



5 Shadowing corrections for the gluon distribution 

To calculate the damping factor using Eq.(44) we need to calculate the SC for the gluon 
distribution. We anticipate that these shadowing correction will be large as 

1. The cross section of the interaction of two gluons with a transverse separation r± 
is bigger than for qq pair and is equal to [9] 



4 r 



2. Each gluon in the parton cascade can interact with the target. Such interactions 
generate the so called "fan" diagrams (see Fig.6), which are essential in the region of 
small X, and correspond to the interaction of a fast GG pair[2] (sec Fig.3). 

In this paper we chose the following strategy of how to calculate the SC to the gluon 
distribution. First, we consider the SC for the fast gluon and discuss the contribution of 
the large distances. After that, we will discuss the contributions of the "fan" diagrams. 

The SC for the fast GG pair has been calculated by Mueller[9], and in our notation 
the gluon distribution has the form { Nc = Nf = 3) 

xG{Q\x) = 1 l^d^r dH^r-d^l^ ^_ ^ _ 

TT^ Jx X' J TT JAr TT r, 

To understand the physical meaning of this equation it is instructive to write down the 
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evolution equation for the gluon density. Indeed, 

d^xG{Q^, x) 



dln{l/x)dlnQ 



2 



(55) 



TT '^^ + (i?'g2)(fc-i) 1^ 2 J ■ 

The first term corresponds to the usual GLAP equations, while the second one takes into 
account the SC. It should be stressed that the term with A; = 1 (if treated as an equation as 

in Ref.[9]), is the same term that appears in the nonlinear GLR equation which sums the 
"fan" diagrams [2]. This term has been calculated using quite a different technique [25, 26]. 
The coefficient in front of the other terms, reflects the fact that all correlations between 
the gluons have been neglected, despite the fact that the gluons are uniformly distributed 
in the target disc with a radius — B' . 

Mueller's formula is not a nonlinear equation, it is the analogue of the Glauber formula 

for the interaction with a nucleus, which gives us the possibility to calculate the shadow- 
ing corrections using the solution of the GLAP evolution equations. Hence, this formula 
should be used as an input to obtain the complete effect of the SC, for the more compli- 
cated evolution equations, such as GLR[2], or the generalized evolution equations [2 7]. 

Galculating the sequent iterations of Mueller's formula gives us a practical way to 
estimate the value of the SG from more complicated Feynman diagrams, that have to be 
taken into account (diagrams such as in Fig.3). Certainly, the iteration of MuUer's formula 
is not the most efficient way to calculate the SC correction in the region of extremely small 
X, but, it gives a systematic way to estimate the value of the SC due to gluon screening 
and it could be a reliable procedure in the HERA kinematic region where x is not too 
small. 

Making use of the exphcit form of S{h±) (see Eq.(29)), we take the integral over 6^ 
and obtain 

xG{Q\x) = l/'^r^^{C + In KG{x',r'l) + ,E,{KG{x'yi) } , (56) 

TT^ Jx X' J-L TT rf 

where C is the Euler constant and Ei is the exponential integral (see Ref.[19] 5.1.11) and 

-oix'yi) = ^-^x'G^'^'{^,,x')). (57) 

In our numerical estimates which will be presented in the next section we use this formula 
to calculate the SC due to gluon screening. We also make the second iteration of Eq.(56) 
substituting the result of the first iteration back into Eq.(56) in place of xG^^^^ . In this 
way we take into account the first contribution of the "fan" diagrams (see Fig. 7). We will 
discuss the result of these calculation in the next section. 
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6 Numerics and comparison with experiment 



6.1 Vector meson production in DIS without SC 

We show in Fig. 8a the dependence of (7(7*p Vp) on x with no SC for some typical 
values of (5^. 

As we have seen in section 2, the pQCD matrix element for DIS production of vector 
mesons is given as a convolution of {Q^ ,r^j z), '^^{rj_,z) and cr{rj_,q±^) where we 
are able to calculate the above when both 7* and V are longitudinally polarized. The 
specifications of ^'^* , and a used in our calculations are as follow: 

1. We use the commonly utilized form for , e.g. our Eq.(lO). 

2. For we use the two forms which were also used in Ref.[7]. The asymptotic form[28] 
denoted AS and the form suggested by Chernyak and Zhitnitsky[29] denoted CZ. We note 
that Braun and Halperin[30] have cast some doubts on the CZ form and have suggested a 
modified wave function not that different from the AS form. The difference between our 
results with the AS and the CZ inputs, as seen in Fig.Sa, defines a theoretical margin of 
error due to the ambiguity in our knowledge of ^f^. This should be compared later on 
with the margin of difference introduced due to the SC. 

3. The qq cross section is given in the LLA utilizing xG'-^^^^, e.g. our Eq.(15), with 
the GRV parametcrization[31]. In our opinion GRV is suitable for our calculations as 
they chose a small scale and as such their result for small x are close to the LLA, enabling 
a reasonably consistent calculation. 

A semiquantitative calculation of (j{j*p pp) has been presented by Brodsky ct.al.[7] 
and is compared to the low energy NMC data[32] a,t x — 0.06 and = 10 GeV'^. We note 
a difference between the approximation used in Ref.[7] and ours. When calculating the 
cross section shown in Fig.Sa we take into account a z dependence of o'{rj_, g_L^) introduced 
after the rj_ integration. Brodsky et.al. neglect this z dependence and as a result their 
calculation with the AS or CZ wave functions have a fixed ratio of (|)^. The equivalent 
ratio in our calculation depends on x and as can be seen in Fig.Sa. Note that in 
our integration we have a cutoff = Q^z > Qq^ = OAGeV"^. The relatively strong 
dependence of the above ratio on x and shows that the LLA used to obtain Eq.(19) 
does not work too well in the kinematical domain under investigation. We shall return 
to this point in section 7. Regardless of the above, the calculation of Ref.[7] and ours at 
X — 0.06 and Q"^ — lOGeV'^ are mutually compatible and are comparable to the NMC 
data[32]. 
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6.2 Vector meson production in DIS with SC 

Fig 8a should be compared with Fig 8b where we present our results for o-(7*j5 pp) 
with SC as defined in section 3. Inspection of Fig.8 shows that the reduction of the cross 
sections due to SC is bigger than the margin of theoretical error shown by the difference 
between the AS and CZ calculations of Fig. 8a. As such SC are demonstrated to be an 
important ingredient of the pQCD calculation and should not be ignored. The comparison 
between the screened and non screened calculations of o"(7*p — > pp) are further shown in 
Fig. 9 as a function of for cm. energies of = 50, 100 and 150 GeV. As expected, 
the difference between the screened and non screened results gets smaller when is 
increased and becomes bigger when is increased. 

Recently, some ZEUS data on (7(7*^ pp) became available [33], and when combined 
with the low energy NMC data[32] provide a first qualitative examination of the pQCD 
results. 

1. The ZEUS differential cross section slope is found to be 5 = 5.1 + 1.2 — O.Oil.OGeV^"^. 
This B value is about a half of the slope found in elastic p° photoproduction and is in 
excellent agreement with our input[13,21]. 

2. ZEUS reports a dependence of (7(7*]? — * pp) where n = 4.2 ± 0.8 + 1.4 — 0.5. 
This should be compared with the (asymptotic) prediction of Brodsky et.al.[7] of n = 6. 
In our calculations n is a function of x and approaching n = 6 when the SC are dimin- 
ished. Note that n — > 2 when In^ is big enough, actually well above the presently available 
HERA kinematics. In the low Q'^ domain investigated by ZEUS, 7 < Q"^ < 25GeV'^, our 
calculations result in n = 5.3 for AS and n = 4.8 for CZ. Clearly the present experimental 
errors are too big for a definite conclusion. Note, also, that the ZEUS results are obtained 
for a = aL + (Jt with R = ^ = 1.5 + 2.8 — 0.6, whereas the pQCD analysis applies to ctl 
only. 

3. The dependence of the ZEUS-NMC cross sections on x is compared with the pQCD 
results in Fig 10 for x < 0.01. The unscreened results calculated by ZEUS according 
to Ref. [7] are given by the shaded area, whereas the screened results are given by the 
two solid lines denoted CZ and AS. At this time the experimental errors are sufficiently 
large so as to accommodate both the screened and unscreened options. However, it is 
encouraging to note that a luminosity increase by a factor of 2-3 will enable a data based 
discrimination between the theoretical models. Note that the theoretical predictions on 
Fig. 10 are obtained after a multiplication by R = 1.5, to account for the observed ar as 
well. 
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Fig. 11 shows the contours of the damping factor D'^ (defined in Eq.(l) and evaluated 
in Eqs. (44,45)) as a function of ^/s and Q"^ for 7*p — > J/iljp. Note that for heavy quarks 
the dependence on is ehminated in the nonrelativistic hmit. Indeed, the difference 
between the AS and CZ resuhs was small enough to be neglected. Experimental data 
on this reaction in the HERA energy range is now available [34] with < AGeV'^. This 
data together with lower energy data is shown in Fig. 12 together with some theoretical 
predictions. The shaded area corresponds to Ryskin's pQCD calculations [5] and our cor- 
rections to Ryskin's predictions are presented by the two solid lines. Note that ours is an 
approximate assessment due to the following uncertainties: 

1. Ryskin uses the BFKL kernel whereas we use a GLAP kernel. 

2. The shaded area in Fig. 11 corresponds to the ZEUS experimental values for xG{x, Q^) 
inserted into Ryskin's calculations. Our calculation of is based on the GRV parame- 
terization [31]. 

3. Ryskin uses a non relativistic quark model wave function for the vector meson, whereas 

we assume the relativistic wave functions which were also used in Ref[7]. This difference 
results in a normalization ambiguity between the two approximations which has not as 
yet been resolved. 

Regardless of the above, an inspection of Figs. 11 and 12 shows that if a pQCD calcu- 
lation, with no SC, yields an energy dependence of s^^ for (T{'~f*p J/ipp), the impo- 
sition of SC will reduce A by approximately 0.08 for the range of Ref.[34]. For real 
photoproduction, = 0, the reduction in A will be 0.11. 



6.3 The damping factor with SC for the gluon distribution 

Fig. 13 shows the damping factor for J/ip production calculated using Eq.(56) for the gluon 
distribution and Eq.(9) for the vector production amplitude. As has been discussed, we 
anticipate a bigger influence of the small distances in the gluon damping than we have 
observed in the the quark sector with SC. The damping factor in the gluon sector is 
calculated to be 

^2 ^ xG{a\z = |),a;)[^^-(^^)l +xG^^^(gg,x) 

where Ql — 0.4 GeV'^. This should be compared with Eq.(44) in the quark sector. In 
Eq.(56) we have integrated over from ^ to Note that the denominator of Eq.(56)is 
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the Born Approximation of the numerator. Since we have seen that the LLA for the GLAP 

evolution equation docs not work well in the kinematic region of interest, we also use this 
approximation for the denominator in our calculation as the SC is known only in LLA. We 
also demand that the initial condition for the gluon distribution is the same for SC and 
without them. Our reasoning is that we look on the initial condition as if it has been taken 
from the experiment. Therefore, the result of the calculation can be read in the following 
way: the suppression of J/'^ production is given by due to the screening at distances 
of 4j > > One can extract from Figs. 12 and 13 that the SC corrections in the 
gluon sector, are as large as in the quark sector. This fact could be anticipated because 
the Kg is | times bigger than Hq. Additional integration over x in Eq.(56) suppresses 
the gluon screening, but not to an extent that will make it smaller than SC in the quark 
sector. 

The ratio R — ^^^Tr"'^'! iteration) sjjQ-^ jj^ Fig. 14. R is the quantitative measure 

Dq( first iteration ) o 

of the correction to the damping factor due to the second iteration of Eq.(56) which takes 
into account the first correction from the "fan" diagrams (see Fig. 7). One can see that 
this correction is rather big and we plan to study them in a more comprehensive way in 
a future publication. 



7 Conclusions 

1. The formulae for shadowing corrections for vector meson production in DIS have been 
obtained whihin the framework of the CLAP evolution equations in the low x region. It 

is shown that the rescattering of the quark is concentrated at small distances ( r_|_ oc ^ ) 
and can be treated theoretically on the basis of pQCD. On the other hand, the rescatter- 
ing of gluons depends on a wider range of the distances, including the larger ones where 
(r± > This causes a large uncertainty in the pQCD calculations. We show that 

the gluon shadowing generates a damping factor which is compatible or bigger than the 
damping in the quark sector. 

2. The numerical calculations of the differential cross sections for vector meson production 
has been performed using two models for the hadron wave function at small transverse 
distances. The results of the calculations show that the corrections due to the SC are 
bigger than the non screened uncertainties. There are two sources for these uncertainties: 
(i) the different nonperturbative models for at r_|_ ^0 which, at the moment, are not 
under full theoretical control. And, (ii) the limited accuracy of the LLA for the GLAP 
equation in the available kinematic region. Basically, those uncertainties can be seen 
without a detailed calculation just by changing the scale of Q'^ (let say, from to Q'^/2). 
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However, we think that the exact integration over z in Eq.(9) has a certain meaning, 
since the GRV gluon structure function reaches the value of the anomalous dimension 
7 = 1/2 at a: ~ 10^^. It means that at low values of x we can try to calculate the gluon 
distribution using the BFKL equation which has a well defined scale in . Therefore, 
the difference in the results due to the integration over allows us to estimate the theo- 
retical uncertainties that originate from our poor knowledge of the nonperturbative wave 
function of the produced vector mesons. 

3. The numerical calculation of the SC shows that (i) they should be taken into account 
in the HERA kinematic region, (ii) their value is bigger than the uncertainties related to 
the unknown nonperturbative part of our calculation, and (iii) DIS vector meson produc- 
tion can be used as a laboratory for the investigation of SC. In our opinion, using these 
reactions to extract the nonperturbative wave function of the produced vector mesons is 
all too ambiguous, at least with the present theoretical knowledge of SC. 

4. The calculation of the SC for the gluon distribution has an intrinsic uncertainty, re- 
lated to the large distance contribution. However, we wish to mention that the anomalous 
dimension in the GRV gluon distribution becomes larger than 1 (7 > 1) at x < 10~^. It 
means that at such small x, all uncertainties from the large distance behaviour become 
rather small, especially at < 2GeV^ ( which corresponds to < lOGeV^ for the 
case of p production and Q'^ < AGeV"^ for J/^ production). We demonstrated that the 
poor accuracy of our calculation in this kinematic region is due to the fact that the LLA 
for the CLAP equation is seriously ill in the HERA region. Much more work is needed 
to provide a reliable calculation of the shadowing corrections for the gluon distribution. 
We established the fact that the Mueller (Glauber - like) formula for the gluon SC cannot 
provide a reliable result as the second iteration of this formula gives a rather big correc- 
tion. We propose to use this formula as an initial condition to solve the new evolution 
equations in the region of small x and we will publish a more detail study of this problem 
in the future. 

5) The relationship between ^^q/Jq2'^^ and the screened cross section for virtual ph- 
toproduction of vector mesons has been established. The relation includes all shadowing 
corrections and in our oppinion is suitable to extract the information on the non pcrtur- 
bative light cone qq vector meson wave function. To do this, we need to know F2'^'\Q'^, x) 
over a wide region of x and Q^. With the appearance of the new HERA data we hope, 
that this will be possible in the very near future. 
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Figure Captions 



Fig.l: Vector meson production in DIS without SC. 

Fig. 2: Extra gluon emission for vector meson production in DIS without SC. 
Fig.3: SC for vector meson production in DIS. 

Fig.4: The time structure of the SC in the inelastic cut of the diagram of Fig.3. 

Fig. 5: The behaviour of the damping factor versus x — Kq for J/-^ diffractive pro- 
duction. 

Fig.6: "Fan" diagrams with triple "ladder" vertices 7. 

Fig. 7: The first iteration of the "fan" diagrams in the gluon distribution for vector 
meson leptoproduction. 

Fig.8: The x dependence of a{-i*p pp) with Q"^ = 4, 9, 15, 20 GeV^. a) In a pQCD 
calculation with no SC. b) In a pQCD calculation with SC. Sohd lines correspond to a 
CZ wave function input and dashed lines to AS. 

Fig.9: Comparison of screened (solid line) and nonscreened (dashed lines) calculated 
values of a{Yp ^ PP) for ^/s = 50, 100, 150 GeV^. a) With AS input, b) With CZ input. 

Fig.lO: (7(7*p ^ pp) as a function of x at = 8.8 and 16.9 GeV'^ taken from Ref.[33]. 
Shaded area corresponds to a pQCD calculation with ZEUS gluon densities and no SC. 
Solid lines are our AS CZ predictions in a pQCD calculation with SC. 

Fig.ll: A mapping of D'^ as a function of ^/s and for J/t/j electroproduction. 

Fig. 12: A compilation of (7(7^ J/i'p) taken from Ref. [34]. The solid line is a VDM 
prediction. The shaded area is the prediction of Ryskin[5] with ZEUS gluon densities. 
The dashed lines are our corrections to Ryskin's calculation. 

Fig. 13: The damping factor for J/il) production due to first iteration screening in 
the gluon distribution for different values of a?. 
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Fig.l4: The ratio R'^ — DQ{second iteration) / Dq{ first iteration) for J/t/: production 
due to screening in the gluon distribution. 
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